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We report on a measurement of splitting in the excitation spectrum of a condensate driven by 
an optical travelling wave. Experimental results are compared to a numerical solution of the Gross 
Pitaevskii equation, and analyzed by a simple two level model and by the more complete band theory, 
treating the driving beams as an optical lattice. In this picture, the splitting is a manifestation of 
the energy gap between neighboring bands that opens on the boundary of the Brillouin zone. 
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Perturbative bulk excitations over the ground state of a 
weakly-interacting Bose Einstein condensate (BEG) have 
been extensively studied using two-photon Bragg transi- 
tions^, 0,0,01 • Their spectrum 5], dynamics 6] and de- 
cay 01 are rather well understood. However, Bragg pro- 
cesses can also strongly excite the condensate, generating 
non-perturbative dynamics, such as Rabi oscillations be- 
tween two (ot more) macroscopically-populated momen- 
tum states 1^ 0, . Rapid oscillations have been found 
to suppress various inhomogeneous broadening mecha- 
nisms, thus increasing the coherence time of the system 

Since the Bragg process may be viewed as the diffrac- 
tion of an atomic wavepacket from a nearly perfect 
optical lattice, band theory may be employed for the 
study of such systems. A BEG in an optical lattice has 
proved to be an excellent system for the investigation 
of effects predicted long ago in solid state physics such 
as Bloch oscillations and Landau-Zener intra-band tun- 
nellin g , fill IT3 . IT3l | . and the Mott-insulator phase tran- 
sition [lj| . Much progress has been made in understand- 
ing the ground state of a BEG in an optical lattice, its 
weak excitations and the effects of interactions [l^ [13 ■ 
However strong excitations involve two beating macro- 
scopically occupied modes, and lead to Rabi like oscilla- 
tions. Thus they cannot be described as a Bogoliubov- 
like excitation over a single Bloch state. 

In this paper we study the spectrum of a BEG driven 
by a strong resonant Bragg pulse, which leads to coher- 
ent Rabi oscillations between momentum states. We ob- 
serve these oscillations in the time domain and a split- 
ting between the energy bands in the frequency domain. 
The splitting is measured experimentally with Bragg 
spectroscopy, using an additional Bragg beam-pair as 
a weak probe, and found to agree with numerical so- 
lution of the time-dependent Gross-Pitaevskii equation 
(GPE). The splitting is also predicted by a simplified 
two-level dressed-state model. However some features are 
explained only by the more complete band theory. These 
theories do not include inhomogeneous and finite-size ef- 
fects, nor interactions, but are in reasonable agreement 
with our experimental results. 

The energy spectrum of a linear system is related to the 



time dynamics via Fourier transform. Therefore the os- 
cillations in the time domain are translated to splitting in 
the spectrum, and the decay of the time-correlation func- 
tion results in broadening of the spectral peaks 
our experiment the chemical potential, which character- 
izes the non-linearity, is smaller than the Rabi frequency, 
but of the same order of magnitude. Hence this linear 
description is only approximately true for our system. 

Our experimental apparatus is described in 5] . Briefly, 
a nearly pure (- 90%) BEG of 5(±1) x 10'' ®'^Rb atoms 
in the \F,nif) = |2,2) ground state, is formed in a mag- 
netic trap with radial and axial trapping frequencies of 
LUr = 27T X 226 Hz and = 27r x 26.5 Hz, respectively. 
The condensate is driven by a pair of strong counter- 
propagating Bragg beams with wave- vectors k^z and 
—ki^z generating an optical lattice potential along the 
axial direction with a depth characterized by a Rabi fre- 
quency The laser frequency is red detuned 44 GHz 
from the ^^Rb D2 transition in order to avoid sponta- 
neous emission. As shown in 0], in the Rabi regime the 
mean field shift is largely suppressed, hence the frequency 
difference between the driving Bragg beams (in the labo- 
ratory frame) is set to w^; = 27r x 15 kHz, the free-particle 
resonance. 

After driving pulses of varying duration, the magnetic 
trap is rapidly turned off. A typical absorption image 
after 38 msec time of flight expansion is shown in Fig. 

The wavepackets with wavenumber — 2/cl and k — 
are clearly separated. In this image, ildt ~ Stt, hence 
most atoms are in the —2kL wavepacket. The number of 
atoms in each wavepacket is measured by integration over 
the density in an area determined by a gaussian fit to the 
absorption image. We thus define iV-2fci and Nq as the 
number of atoms in the wavepackets with wavenumbers 
— 2fci and respectively. 

Figure 121 shows the time dynamics of the BEG exposed 
to the driving Bragg beams. The solid squares are data 
points, measuring the fraction of atoms in — 2fci cloud 
as a function of the duration of the driving pulse. The 
BEG undergoes oscillations at a frequency of 27r x 5.3 ± 
0.1 kHz as found from a sinusoidal fit to data. The hollow 
points are numerical solutions of the GPE including finite 
time, and spatial inhomogeneity with Qd as the only fit 
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FIG. 2: Rabi oscillations between modes with wavenumbers 
and -'2kL (in the laboratory frame of reference). The solid 
line is a sinusoidal fit to the data which gives an oscillation 
frequency of 2n x 5.34 kHz. The open circles are the GPE 
simulation with fid = 27r x 5.5 kHz whose sinusoidal fit gives 
the same oscillation frequency. The oscillation frequency is 
less than fid due to interactions, which in effect lower the 
depth of the optical lattice. 



FIG. 1: Time of fiight images of BEG undergoing strong 
Rabi oscillations, (a) BEG evolved in the presence of driving 
Bragg beams alone. The pulse of duration 320 jisec is just 
over 37r/r2d, leaving almost all population in —2kL momen- 
tum state, (b) Both driving Bragg beams and probe Bragg 
beams are on for 600 fisec. Driving pulse is close to 6n/Q,d 
leaving the state — 2A;l almost unoccupied. 



parameter [l^]. The Rabi frequency which yields the 
correct oscillation frequency is found to be fid — 2ttx 5.5 
kHz. The oscillation frequency is slightly less than ild 
due to interactions. This is because the higher density in 
the valleys of the optical potential creates a mean-field 
lattice potential opposing the optical one 0| . 

Next, we probe the spectrum of the oscillating BEC, 
with another pair of Bragg beams (turned on simultane- 
ously with the driving Bragg beams), generating a weak 
lattice with a Rabi frequency flp. The probe beams, also 
counter-propagating along the axial direction, are lin- 
early polarized perpendicular to the driving beams. The 
frequency difference between probe beams LUp is always 
kept such that the momentum wavepacket is diffracted 
by the probe lattice to a 2kL momentum wavepacket 
whereas for the — 2fci wavepacket, the Bragg-frequency 
of the probe beams is always doppler-shifted away from 
resonance. For the same reason, the driving beams are 
Doppler-shifted from the state. A characteristic time 
of flight image is shown in Fig.^. Here in addition to the 
clouds in states — 2fci and 0, there are also atoms with 
wavenumber 2kL. The measured relative population of 
this momentum state N2kL/{^-2kL + + -^2fci) as a 
function of the probe frequency difference Up forms the 
probe Bragg spectrum, which is plotted in Fig. |31 (black 
squares). The duration of the beams is set so the driving 



pulse would be close to an even multiple of tt, in order to 
minimize collisions in the time of flight expansion. 

The probe spectrum, without the driving beams {fid — 
0) is similar to the well known Bogoliubov spectrum 
(Fig.|21i). The probe is weak in the sense that ftp ~ f^d/S, 
but since the probe pulse is greater than 7r/2 we observe 
a suppression in the mean field shift Q . The probe spec- 
trum in the presence of the driving beams (Fig. ^p) 
shows a clear splitting in the spectrum of the oscillating 
condensate. The dotted line is a double Lorentzian fit 
to the experimental data, with a distance between the 
peaks of AE/h = 4.7 zL 0.2 kHz, comparable but larger 
than the chemical potential ii/h = 1.6 kHz. The solid 
line is the full GPE simulation of the experiment, which 
gives a splitting of AE/h = 5.5 kHz. We note this is 
not the frequency of oscillations in Fig.|21 but rather the 
Rabi frequency of the lattice itself. The mismatch be- 
tween experimental and simulation results, may be due 
to a drift in the intensity of the laser beams between the 
oscillations measurement and the splitting measurement. 

The width of the spectrum without driving (Fig. |3^) 
fits the GPE simulation, and is mostly due to finite time 
of the experiment. However the experimental peaks in 
Fig. are substantially broader than predicted by the 
GPE. This broadening is robust and is observed for dif- 
ferent values of Qd, indicating broadening mechanisms 
which are beyond mean field. One such mechanism for 
a homogeneous condensate is due to uncertainty in a 
dressed state decomposition of the initial Fock state, that 
leads to a spectral broadening 19J. For our experimental 
parameters this broadening is much less than the widths 
of the experimental peaks which are 3 ±0.4 and 3.7 ±0.5 
kHz. We do not fully understand the broadening. The 
overall response of the system to the probe is smaller in 
Fig. |3Jd (with driving beams) than that in Fig. |2K (with- 
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FIG. 3: Excitation spectra: Fraction of atoms in 2kL state 
as a function of the probe frequency difference, (a) No lattice 
{Q.d = 0) and (b) fid = 5.5 kHz. In (a) The solid line is a 
GPE simulation with Q.p as a single free parameter. In (b) 
the solid line is a GPE simulation with f2p the same as in (a) , 
and f2d determined by Fig. |21 The dotted line is a Lorentzian 
fit to the data giving a splitting of AE = h x 4.7(±0.2) kHz. 
The dashed line is the noninteracting band model discussed 
in the text. 

out driving beams). This is because the time-averaged 
population of the state fc = in the presence of the lattice 
is only half of that of the ground state of the BEC. 

The simplest model that explains the splitting in the 
spectrum is a noninteracting model of a closed system 
with two momentum states: and 2hkL, coupled to a 
classical light potential. Choosing our frame of reference 
as moving with the optical lattice of the Bragg driving 
beams at velocity v = —k^lM along z yields a time- 
independent Hamiltonian, 



Ha = -77T7^ + ^^dCOs{2kLz). (1) 

The eigenstates of such a Hamiltonian are dressed 
states |20|. The energy difference between the two 
eigenenenergies is firf. This energy difference is mani- 
fested both in the splitting of the spectrum, and in the 
Rabi oscillations which are the coherent beating between 
the two eigenstates. Since the initial k — Q state is an 
equal super-position of the two dressed states, the two 
peaks in the spectrum are expected to have the same 
hight. 

The system is however not a closed two-state system, 
as there are nonresonant Rabi oscillations with other mo- 



mentum states. We therefore calculate the eigenstates 
of the non-interacting Hamiltonian ^ without the re- 
striction to two momentum states. According to Bloch's 
theorem, due to the periodicity of the optical lattice, a 
state \k) with wavenumber k is only coupled to states 
\k + 2mki,), where m is an integer, and 2kL is the basis 
of the reciprocal lattice. The Hamiltonian is diagonalized 
by the Bloch states 

|n, <7> ^ an^q(k)\q + 2mkL) (2) 

m 

where q is the quasi-momentum and n is the band index 
0. In the moving frame of reference, the stationary BEC 
has momentum hk^ and is situated on the Brillouin zone 
boundary. The initial kinetic energy of the condensate is 
in the energy gap, and by suddenly turning on the light 
potential, the initial state jfei) is projected onto the Bloch 
states |fci) — X^n '^"(^i)l'^)' where we have omitted the 
quantum number q which is understood to he q = kL- 

In the weak lattice limit, hfld <ti Er where Er — 
f{^k\/2M is the single photon recoil energy, we recover 

the two state result jfe^) — l/v^(|l) + |2))i and a split- 
ting of For our parameters the band structure is 
obtained by numerical diagonalization and presented in 
Fig. 01 The predicted splitting between the two lowest 
bands is indeed within 1% of fid, but contrary to the weak 
lattice approximation, \a^{±ki)\ = \a^^{±kL)\ ~ 10% of 
|a3(±3fcL)| = |a4(±3A:L)|. This indeed leads to negligi- 
ble occupations, but as we show shortly, to observable 
interference effects. 

With the probe beams, the Hamiltonian governing the 
model homogeneous system in the driving lattice frame 
is H = Hq + Hp with 

Hp = hQp cos{2kLZ - ujpt) = Fpe"'"^"* -t- c.c. (3) 

with Vp — MlpC^^^'^^ 12. In the rotating wave approxi- 
mation Hp is a momentum shift operator coupling state 
|fc) to |fc -I- 2/cl). Since Q.p « fi^, it can be treated by 
the use of first-order time-dependant perturbation the- 
ory, giving probability amplitudes for transitions between 
bands. The momentum transferred to the condensate by 
the redistribution of photons in the probe beams is a mul- 
tiple of the reciprocal lattice vector , hence also with 
the probe beams the initial quasi-momentum q = kL is 
conserved. 

The dashed line in Fig. is the full first order per- 
turbation theory expectation of the population of state 
|3fci) in the moving frame of reference, which is equal to 
2kL in the laboratory frame of reference. Aside from a 
mean- field shift, this simple model seems to fit the CPE 
simulation quite well. 

Since the splitting between band 3 and 4 is negligible, 
the splitting in the spectrum of the probe in this non- 
interacting model is essentially the splitting between the 
two lowest bands (see Fig. 0^). 

There is a clear asymmetry in the peak heights which is 
in good agreement with the CPE simulation. This asym- 
metry is due to quantum interferences between pathways 
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FIG. 4: Noninteracting band structure for our experimen- 
tal parameters, (a) The energy of the four lowest bands as a 
function of quasi-momentum (solid lines) . The energy gap on 
the boundary of the Brillouin zone is 0.99 hQ,d, while the gap 
between the third and fourth band is two orders of magnitude 
smaller. The dashed line is the free particle energy, (b) The 
amplitudes of the different momentum states k — kL + 2mkL 
in each Bloch state on the boundary of the Brillouin zone. 
Black, dark gray, light gray and white are bands 1,2,3 and 4 
respectively. The arrows represent the coupling by the probe. 
Note the contamination of each Bloch state by the "unintu- 
itive" momentum states. 

of the probe excitation. The probe couples each of the 
lower bands (1 and 2) to each of the higher bands (3 
and 4) by two paths. In the weak lattice approxima- 
tion there is only the path \ki,) |3fci), however due 
to the nonvanishing az{kL) and aA{kL), there is also a 
path I — /cl) ^ I^l). These paths are marked by arrows 
connecting the relevant amplitudes in Fig 4(b). The in- 
terference between both paths is constructive for tran- 
sitions from band 2 and destructive for transitions from 
band 1. Therefore the lower energy peak in the split- 
ting, corresponding to transitions from band 2 should be 
larger than the higher energy peak which corresponds to 
transitions from band 1. 

In the experiment, the frequency difference between 
the driving Bragg beams is equal to the free particle res- 
onance fvjJd = ^Er. For other values of fvaJd, the Rabi 
oscillations are not full. As the detuning from resonance 
grows, the amplitude of oscillations is reduced and the 
frequency of oscillations grows. In the spectral domain 
the frequency difference between the two peaks of Fig. 

grow with the detuning from resonance. The ampli- 
tude of one peak grows while it's position in frequency 
shifts towards the excitation energy of a free Bogoliubov 
quasiparticle. The other peak shrinks in size as the de- 



tuning of ujd from resonance grows. In Fig. [Slwe plot the 
values of ujp at the peak of the spectrum versus uJd, as ob- 
tained by numerically solving the GPE for Vld = 5.5 
kHz. We compare this to the energy splitting of a non- 
interacting two-level system, calculated for the same pa- 
rameters, and shifted up by 1.3 kHz. As in the resonant 
case, the splitting is in good agreement with the GPE, 
except for a constant shift due to interactions. 
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FIG. 5: The locations of the peaks in the response to the 
probe, as a function of the driving beams' frequency differ- 
ence LOd- (circles) - numerical solution of the GPE including 
finite size, finite time and mean field effects, (lines) - the en- 
ergy levels of a two-level system shifted up by 1.3 kHz. The 
probe wavenumbers are taken as ±fc_L/3 to distinguish probe 
excitations from non-resonant driving excitations. 

In conclusion, we measure a splitting in the Bragg ex- 
citation spectrum of a BEG undergoing Rabi oscillations 
between two momentum states. Experimental data fits 
well with numerical GPE simulations. The main features 
are captured by a simple noninteracting model, whereas 
the main contribution of interactions is found to be the 
mean field shift of the probe spectrum, captured by the 
GPE, and an increase in the width which is not cap- 
tured by the GPE. The splitting in the spectrum results 
in a splitting in the energy of the atoms undergoing inco- 
herent collisions. Experimentally, a divergence from the 
known s-wave scattering halo was measured 21] for a 
rapidly oscillating BEG. There too, a pronounced asym- 
metry was found between the inner coUisional shell and 
the outer one. 

This work was supported in part by the Israel Min- 
istry of Science, the DIP foundation, the Israel Science 
Foundation and by the Minerva foundation. 



[1] F. ZambeUi, L. Pitaevskii, D. M. Stamper-Kurn, and [2] P. B. Blakie and R. J. Ballagh, J. Phys. B 33, 3961 
S. Stringari, Phys. Rev. A 61, 063608 (2000). 



5 



(2000). 

[3] J. Stenger, S. Inouye, A. Chikkatur, D. Stamper-Kurn, 

D. Pritchard, and W. Ketterle, Phys. Rev. Lett. 82, 4569 
(1999). 

[4] D. Stamper-Kurn, A. Chikkatur, A. Gorlitz, S. Inouye, 

S. Gupta, D. Pritchard, and W. Ketterle, Phys. Rev. 

Lett. 83, 2876 (1999). 
[5] J. Steinhauer, R. Ozeri, N. Katz, and N. Davidson, Phys. 

Rev. Lett. 88, 120407 (2002). 
[6] R. Ozeri, J. Steinhauer, N. Katz, and N. Davidson, Phys. 

Rev. Lett. 88, 220401 (2002). 
[7] N. Katz, J. Steinhauer, R. Ozeri, and N. Davidson, Phys. 

Rev. Lett. 89, 220401 (2002). 
[8] J. Hecker Denschlag, J. E. Simsarian, H. Haffner, 

C. McKenzie, A. Browaeys, D. Cho, K. Helmerson, S. L. 

Rolston, and W. D. PhilUps, J. Phys. B 35, 3095 (2002). 
[9] N. Katz, R. Ozeri, E. Rowen, E. Gershnabel, and 

N. Davidson, Phys. Rev. A 70, 033615 (2004). 
[10] A. S. Mellish, G. Duffy, C. McKenzie, R. Geursen, and 

A. C. Wilson, Phys. Rev. A 68, 051601 (2003). 
[11] O. Morsch, J. Miiller, M. Cristiani, D. Ciampini, and 

E. Arimondo, Phys. Rev. Lett. 87, 140402 (2001). 

[12] F. CataUotti, S. Burger, C. Fort, P. Maddaloni, F. Mi- 
nardi, A. Trombettoni, A. Smerzi, and M. Inguscio, Sci- 
ence 293, 843 (2001). 

[13] B. Anderson and M. Kasevich, Science 282, 1686 (1998). 

[14] M. Greiner, O. Mandel, T. Esshnger, T. Hansch, and 
I. Bloch, Nature 415, 39 (2002). 



[15] M. Kraemer, C. Menotti, L. Pitaevskii, and S. Stringari, 
Eur. Phys. J. D 27, 247 (2003). 

[16] D.-I. Choi and Q. Niu, Phys. Rev. Lett. 82, 2022 (1999). 

[17] D. Tannor, Introduction to Quantum Mechanics: A 
Time- dependent Perspective (University Science Books, 
Sausalito, 2005). 

[18] We numerically solve the GPE for the order param- 
eter of the condensate as described in [53 ihdtip = 
{-fi^V^/2M -I- y -I- gW'^}il^, with the time depen- 
dant external potential typically of the form V{r, t) = 
M/2(Li;^r^ +Lo'^z'^) + hQdCos{2kLZ — LOdt) + hflpCos{2kLZ — 
ujpt). All other parameters for the simulation are as in 
the experimental system, described above. The ground- 
state is found for t < hy imaginary time evolution. We 
exploit the cylindrical symmetry to evolve the wavefunc- 
tion on a two dimensional grid A''^ x A'r(4096 x 32), using 
the Crank-Nicholson differencing method with alternat- 
ing direction implicit algorithm. 

[19] E. Rowen, R. Ozeri, N. Katz, R. Pugatch, and N. David- 
son (2005), cond-mat/0506318. 

[20] C. Cohen- Tannoudji, J. Dupont-Roc, and G. Grynberg, 
Atom-Photon interactions (John Wiley & Sons Inc., 
1992). 

[21] N. Katz, E. Rowen, R. Ozeri, and N. Davidson (2005), 

cond-mat/0505762. 
[22] W. Bao, D. Jaksch, and P. A. Markowich, J. Comput. 

Phys. 187, 318 (2003). 



